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We demonstrate theoretically the feasibility of selective self-excitation of higher-mode flexural 
vibrations of graphene nano-ribbons and carbon nanotubes by the means of magnetomotive insta- 
bility. Apart from the mechanical resonator, the device consists only of a constant voltage source, 
an inductor, a capacitor, a gate electrode and a constant magnetic field. Numerical simluations were 
performed on both graphene and carbon nanotubes displaying an overall similar behaviour, but with 
some differences arising mainly due to the non-linear mechanical bending forces. The advantages 
and disadvanatges of both materials are discussed. 



I. INTRODUCTION 

The interplay between the mechanical and the elec- 
tronic degrees of freedom of nano sized devices continue 
to be an active research area on both the experimental 
and theoretical levels, largely due to the discovery of 
new high performing materials such as graphene and 
carbon nanotubes'. Interest lies not only in the possible 
technological applications of such materials but also on 
a more fundamental level, exemplified for example in the 
numerous schemes for mechanical ground state cooling 
that have been proposed in recent years. 

In two previous papers 3,4 it was explored theoreti- 
cally the electro-mechanics of a carbon nanotube that 
was incorporated as a displacement sensitive resistor in 
a simple electronic circuit. The two device geometries 
considered were an RC-circuit under constant current 
bias and an RLC-circuit under constant voltage bias. 
It was shown that under the influence of a constant 
magnetic field the electronic circuit would either pump 
energy into the mechanical system, possibly leading 
to an instability and selfexcitation of vibrations, or 
the opposite, to drain the mechanical oscillation of 
energy leading to an effective cooling. Remarkably, in 
most cases, by simply reversing the direction of the 
magnetic field one could switch between self-oscillation 
and cooling. On the most basic level, the difference in 
performance between the two different device geometries 
consisted in the fact that in the voltage-bias regime the 
driving voltage could be lowered, but at the expense 
of having to include a sufficiently large inductance in 
the system. On a more abstract level, an interesting 
qualitative difference showed up due to the different 
dimensions of the corresponding dynamical systems. In 
an RC-circuit there is essentially only one electronic time 
scale involved, namely the RC-frequency a;_R ~ 1/{RC). 
In an RLC-circuit, howev er, we also have the LC- 
frequency ujl = \/l/{LC). The consequence of this 
extra time constant turned out to be a more complicated 
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dependence of the overall performance on the relative 
magnitudes of the three frequencies wq, ujr and wl, 
where luq is the mechanical resonance frequency. Most 
importantly, not only the direction of the magnetic field 
determined whether there would be pumping or cooling 
in the system but instead the direction of the magnetic 
field in combination with the relative magnitudes of 
iJL and Wo- To be more precise: Suppose that you 
first fix the direction of the magnetic field so that it 
causes either pumping or cooling. If you then adjust 
the LC-frequency so that ujl — ujq changes sign then the 
effect is invariably reversed, that is, pumping switches to 
cooling and vice versa. In this paper we will show that 
this feauturc, more than beeing a mere curiosity, could 
actually become useful. 

In the previous work only the fundamental bending 
mode was taken into account. However, by considering 
also the higher harmonics a graphene sheet or a carbon 
nanotube could be viewed as an infinite set of mechanical 
oscillators with different frequencies. It is also likely 
that the amplitude of each of these modes would have 
an effect on the resistance. If we again consider the 
voltage-biased device, depicted in Figure (1(a)), and 
adjust the LC-frequency so that luq < lul < ^2, where 
u]2 is the frequency of the second harmonic, by adjusting 
the direction of the magnetic field so that it pushes the 
nanotube or graphene towards the gate electrode we 
could achieve a situation where the second harmonic 
is self-excited while at the same time the fundamental 
mode is kept silent (or even cooled down). This could of 
course be generalized, if you instead want to selectively 
self-excite the 2n:th mode you adjust the LC-frequency 
so that W2n-2 < < ^2n (The reason why only every 
second mode can come into consideration will become 
clear later). For previous results on selective excitation 
of higher vibrational modes see for example Refs 5,6. 
The question is now how far you can actually push this 
scheme in practice. In principle there is no upper limit to 
the higher mode frequencies, extending for example into 
the long sought after THz regime. There are however a 
number of limiting factors. First of all the effect of the 
Lorentz force declines as we go up into higher modes. 
Secondly, it is plausible that the sentivity in resistance to 
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higher mode deflections dechne at more or less the same 
rate. The latter fact also has the consequence that the 
amplitude of the current oscillations generated by the 
mechanical oscillation is reduced accordingly. Another 
complication arises from the non-linear forces which 
become increasingly important at large amplitudes and 
for higher modes. 

A conclusive answer to how high frequency you can 
reach is of course hard to establish, especially since the 
quality of the materials is continually enhanced. The aim 
of this paper is instead to demonstrate through numerical 
simulations the validity of the basic conjecture, namely 
the possibility to selectively self-excite higher modes. For 
each mode considered we will try to assess appropriate 
boundaries for certain parameters, especially the mag- 
netic field strength. Because of the generality of the 
scheme there is no reason to exclusively treat cither car- 
bon nanotubes or graphcne nano-ribbons, instead we will 
discuss and compare the advantages and disadvantages of 
both materials. 



II. THE ELECTRONIC CIRCUIT 
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(a) Sketch of the proposed electronic circuit. 




(b) Sketch of the mechanical resonator comprising the resistor. 

FIG. 1: Sketch of the proposed electronic circuit (upper fig- 
ure). The resistor is comprised of either a grapheme sheet 
or a carbon nanotube suspended over a gate electrode (lower 
figure) . 



The electronic scheme is depicted in Figure (1(a)). The 
resistor is comprised of either a carbon nanotube or a 
graphcne nano-ribbon that is suspended over a gate elec- 
trode. The equations governing the electronic subsystem 
reads: 



CV = I -V/R[z{t,x)] 
Ci = Va-V 



(1) 
(2) 



The resistance R is now treated as a functional of 
the total vertical bending shape z{t,x) of the nan- 
otube/graphene where a; is the horizontal coordinate, see 
Figure (1(b)). The selective self-excitation scheme re- 
quires that each mechanical bending mode of the mate- 
rial affects the resistance. The way in which this occurs 
could be outlined as follows. Let the charge per unit 
length on the resistor be given by the expression 



n = VGCGizit,x)), 



(3) 



where Vq is the gate voltage and CG{z{t,x)) is the gate 
capacitance per unit length. (This could be thought of 
as the defining equation for the gate capacitance). For a 
graphcne sheet the following expression for the classical 
capacitance per unit length was derived in Ref 7: 
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(4) 

where w is the width of the sheet and d is the distance 
to the gate electrode. In the limit w ^ d the above 
expression simplifies to that of a plate capacitor: Cq = 
TreoCrw/d, and in the limit w <^ d it becomes Cg = 
2TTe/ \n{4d/w). The latter expression is similar to that 
for a single walled carbon nanotube: 
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(5) 



which holds in the limit r <^ d, where r is the radius 
of the tube. When we come to discuss the dependence 
of the resistance per unit length on the charge carrier 
concentration n there is a difference. For a single walled 
semiconducting carbon nanotube, for sufficiently low car- 
rier concentration then to a first approximation we have 
that^-i" 



whereas for gra 
hold" 



{C('^)}swNT 2 ' 

iphene we expect the following relation to 

{C("-)}GraphGnc ~ (7) 



In our simulations the total resistance is calculated ac- 
cording to the formula 



R[z{t,x)] = / an{CG{z{t,x))))dx, 



(8) 
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hence, each segment is treated as an independent resistor 
whose capacitance only depends on the distance to the 
gate. This is of course only an approximation. However, 
we have no reason to believe that the exact quantitative 
description of the resistance affects the main qualitative 
conclusions of our analysis, and it will therefore not be 
discussed at length here. For further details we refer to 
the references provided. 



III. CARBON NANOTUBES 

We will first carry out the analysis on a carbon 
nanotube, which was the mechanical oscillator first con- 
sidered. The generalization to a graphene nano-ribbon 
is then straightforward. 

If wc model the carbon nanotube as a doubly clamped 
elastic beam of length L and include a geometrical non- 
linearity term and the external Lorentz force Fl = HJ, 
where H is the magnetic field and J the current through 
the nanotube, the equation of motion reads 



ES 



pA 



dz 




HJ. 



(9) 



Here E is the Young's modulus, p the mass density, A the 
cross sectional area, S the area moment of inertia, 7 the 
damping coefficient and L the length of the nanotube. 
If we let r denote the radius of the nanotube we have 
that A = ■nr^ and 5' = nr'^ /A. For previous work on the 
nonlinear dynamics of carbon nanotubes see for example 
Ref 12. We chose to treat the problem numerically using 
a Galerkin reduced-order model, which we truncate at 
the 6th overtone. Hence 



z{t,x) = w„(t)(/)„(x/L) 



(10) 



ri=0 



Using the notation x 
by the expression 



x/ L, the mode shapes are given 



<j)n{x) =C„{(sin(fc„) - sinh(fc,i))(cos(A:„i) - cosh(fc„a;)) 
— (cos(/s„) — cosh(fc„))(sin(fc„i;) — sinh(fc„i:)}, 

where C„ are normalization constants chosen so that 
4>n{x)^dx = 1 and the constants kn satisfy the equa- 
tion cos (A;„) cosh (A;„) = 1. The corresponding vibrational 
frequencies are given by 



(11) 



Setting the timescale to t = Lo^t and projecting equa- 
tion (9) onto each different mode we obtain the set of 




equations 

Unir) -|-7U„(t) -I- (w„/wo)^Un(T) = Un(T)x 

where 

Xn^ {(j}n{x)(j)'l,{x))dx 
Jo 



K = k^ES/L^ 

an = / (j)n{x)dx 



LHJ 

' K\ 
(12) 



(13) 

(14) 

(15) 
(16) 



The constant Q = I/7 will be referred to as the quality 
factor of the nanotube. We have assumed that the deflec- 
tion is small enough that the Lorentz force can be con- 
sidered uniform across the nanotube. In particular this 
implies that the external force vanishes for n = 1,3,5... so 
that we only need to include the modes n = 0,2,4,6 into 
our consideration. In Table (I) are listed some of the 
mode dependent parameters. We will now review some 



TABLE I: Mode dependent vibrational frequencies cjn and 
integration constants a„ for the first four bending modes with 
even index of a carbon nanotube. 



n 









1 


0.83 


2 


5.4 


0.36 


4 


13.3 


0.23 


6 


24.8 


0.17 



of the analytical results that will form the basis of rea- 
soning in this paper. The measure of the sensitivity of 
the resistance with respect to displacement was expressed 
by the so called characteristic length, which in order to 
include higher modes can be generalized in the following 
way: 



An — 



R 



dR/dur, 



(17) 



In the case which we consider, when the sensitivity of 
the resistance is achieved through electronic doping con- 
trolled by a gate electrode, it could be argued that the 
relative magnitude of the characteristic lengths should be 
approximately 



(18) 



We now introduce the mode dependent coupling param- 
eter 



HJL 



(19) 
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and the mode dependent succeptibility function"' 

UJnUJR{cjl/ojl - 1) 



I) 



(20) 



Neglecting non-linear terms and other intermode 
crosstalk it can be shown that if the product jSnSn is neg- 
ative the mode becomes unstable approximately when 



\PnSn\ > l/Q- 



(21) 



As mentioned earlier, if for example we assume that 
^271-2 < uiL < uj2n then for a fixed magnetic field it 
is clear that P2n-2S2n-2 and P2nS2n have opposite 
signs. If finSn is negative then it will also be negative 
for higher modes but since the averall succeptibility 
decreases with n it should be possible to find a magnetic 
field strength where the slight heating effect on the 
higher modes does not turn into an instability. What we 
have discussed so far is a more or less straightforward 
generalization of previous considerations. The main 
practical complication that arises here is due to the 
possible interaction between the modes mediated by 
the non-linear force term. On this issue, the radius of 
the nanotube is one of the important parameter that 
we need to take into account. Since the nonlinear-force 
becomes appreciable only when the amplitude reaches 
a magnitude of the order of the radius one might think 
that a large radius is beneficial for our purposes. Indeed, 
a large radius not only increases the range of amplitude 
within the linear regime but also increases the spring 
constant K with the side effect that the thermal noise 
is reduced, which can be seen from the equipartition 
relation l/2kBT = l/2K{u^). But, this comes at the 
cost of a reduced coupling parameter. If we make the 
estimates^ ' E = 1 TPa, p = 1.35g/cm^, L = Ifim, in 
Table (II) we have calculated the key parameters for a 
number of different radii. As we can see, going from 



TABLE II: Effective spring constant, thermal noise and vi- 
brational frequency of a carbon nanotube for three different 
radii. 



r [nm] 


K 


[nm] at T=1K 


iJo [lOVs] 


0.5 


2.45-10"^ 


0.75 


0.15 


1.0 


3.93-10-* 


0.19 


0.30 


5.0 


2.45-10"^ 


0.0075 


1.52 



a radius of 0.5 nm to a radius of 5 nm means that K 
increases by four orders of magnitude, and consequently, 
the critical magnetic field increases accordingly. 

Apart from the excitation of an overtone, one can also 
anticipate a certain cooling of the lower harmonics. In 
case we wish to selfexcite the mode with frequency u)2n 
the choise of uj^ and tun that simultaneously maximizes 



the excitation succeptibility S2n and the cooling succep- 
tibility 5*271-2 can be derived straightforwardly: 



(W2„ 



-^2n-2)/i^2n + W2«-2) 



'-^2n-2^B. + ^2n-2 



(22) 
(23) 



In our simulations we took these as our default values 
for the electronic frequencies. We chose to consider a 
nanotube with radius r — Inm and length L = 1/zm 
situated at a distance 200 nm from the gate and put 
A=lnm. Numerical simulations were performed using a 
Runge-Kutta algorithm for different values of Q, uir, ujl 
and the coupling parameter 



HVgL 

KRoX' 



(24) 



where i?o is the resistance of the unbent nanotube. In 
each time step the normalized resistance was calculated 
through numerical integration using 100 steps. The algo- 
rithm was implemented in both Matlab and C producing 
consistent results. 
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FIG. 2: Saturation amplitudes Ao, A2, A4 and Ae for the 
respective modes ito, U2, 1*4 and ue, as a result of computer 
simulations aimed at selective self-excitation of the second 
harmonic of a carbon nanotube with radius r = Inm and 
length L — 1/im. The coupling parameter /? was in the range 
-6.6 to -7.1 which, assuming a stationary current of 1 ^A, 
corresponds to a magnetic field strength in the interval 2.6- 
2.8 T. The other parameters were Q = 10*, ujl ~ 2.323, 
uh = 4.397. Initial conditions were in all variables. 

In Figure (2) we show the saturation amplitudes Ao, 
A2, A4 and Aq for the respective modes uq, U2, U4 and 
uq, as a result of computer simulations for a coupling 
parameter /? in the range -6.6 to -7.1 which, assuming a 
stationary current of 1 fxA, corresponds to a magnetic 
field strength in the interval 2.6-2.8 T. The quality 
factor was assumed to be Q = 10^ and the electronic 
frequencies where chosen so as to selectively excite the 
second overtone. Figure (3) shows the envelopes of the 
rapid oscillations as a function of time of the zeroth and 
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FIG. 3; Time evolution of the envelopes of the rapid oscilla- 
tions in the fundamental mode Mo (upper figure) and second 
harmonic U2 (lower figure) as the results of a computer simu- 
lation on a carbon nanotube with radius r — Inm and length 
L — Ifim. The dashed lines mark the displacement aver- 
aged over a period of 2n. The parameters were /3 = —0.007, 
Q = fO**, UJL ~ 2.323, iJR, = 4.397. Initial conditions were 
in all variables. 



second mode for one of these simulations. As we can see, 
in the beginning the fundamental mode oscillates with 
the largest amplitude but later on the second harmonic 
experiences a rapid boost which causes the magnitude of 
the timeaveraged displacement in the fundamental mode 
to decrease. This is because the oscillation increases the 
effective resistance leading to a reduced Lorentz force. 
The latter phenomenon was explained and quantified 
in Refs 3,4. The other modes also experience a certain 
excitation but with an amplitude roughly two orders of 
magnitude less than that of the second harmonic. 

There are a number of ways by which the non-linear 
forces can cause the scheme to behave in an erratic fash- 
ion. If the amplitude of some mode is forced beyond a 
certain critical value, non-linear forces causes the system 
either to enter a chaotic state which cannot be tracked 
numerically, (which was the typical scenario for carbon 



nanotubes), or to collapse and start all over again, (which 
was the more common outcome for graphene sheets, to 
be disussed later) . Another effect of the non-linear forces 
is that the static deflection increases the effective spring 
constant so that one might end up with a situation where 
the renormalized frequency of the fundamental mode is 
suddenly larger than lul. Although the latter effect is 
only an issue when wl is relatively close to wq , also when 
attempting to self-excite higher modes the static deflec- 
tion induced by the Lorentz force might cause a signifi- 
cant decrease in the effective electro-mechanical coupling 
so that for quality factors below a certain value there 
might not be any possibility for self-excitation whatso- 
ever. For example, given a quality factor of Q = 10''' we 
were not able to find any range of (3 corresponding to 
a successful self-excitation of the fourth harmonic. How- 
ever, by increasing the quality factor slightly to Q = 2-10'' 
tliis was achieved in a number of computer simulations 
presented in Figure (4). As in the previous case, the sat- 
uration amplitude of the other modes was approximately 
two orders of magnitude less than that of the excited 
mode. Assuming, as before, a stationary current of 1 fiA 
the parameter range in this case corresponds to a mag- 
netic field in the interval 9.0-10.6 T. 

Q = 2-1 0^ CO, = 8.485, co„ = 7.943 
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FIG. 4: Saturation amplitudes Aq, A2, A4 and A^ for the 
respective modes uq, U2, U4 and Me, as a result of computer 
simulations aimed at selective self-excitation of the fourth har- 
monic of a carbon nanotube of radius r — Inm and length 
L = Ifim. The coupling parameter /3 was in the range -23 to 
-27 which, assuming a stationary current of 1 fiA, corresponds 
to a magnetic field strength in the interval 9.0-10.6 T. The 
other parameters were Q — 2 ■ 10*, ljl ~ 8.485, ljr — 4.397. 
Initial conditions were in all variables. 



IV. GRAPHENE 

We will now carry out the same analysis on a graphene 
nanoribbon. The elastic and dynamic properties of 
graphene has been analyzed by several authors^ . Af- 
ter a few simplifications, for example that we only need 
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to take into account the vertical streching, the one- 
dimensional equation of motion of a doubly clamped 
graphene sheet reads 

d'^z{t,x) dz{t,x) d'^z{t,x) 

where p is the area mass density of graphene, Pz{t,x) 
is the pressure in the vertical direction, Tq = (A + 2/j,)d, 
Ti = X/2 + IJ,, A and fj, beeing the Lame-parameters and S 
a parameter that quantifies the initial in-plane streching. 
The important thing to notice is that the nonlinear term 
is not dependent on the width of the sheet, unlike the 
carbon nanotube where the radius entered explicitly 
into the equation. Here, it is instead primarily the 
length of the nanoribbon that determines the critical 
amplitude of motion when non-linear forces become 
dominant. In our simulations we chose the typical 
parameter values A — 15.55 J/m^, /i — 103.89 J/m^, 
S = 0.5%^', which together with the length L = 1/im 
yields a fundamental frequency of the order of IGHz. 
For the sake of comparison with the previous simula- 
tions on carbon nanotubes we chose the width w ~ lOnm. 

As before we assume a vertical displacement of the 
form: 

6 

z{t,x) = X^Un{T)(j)n{x) (26) 
n=0 

where the mode shapes are now given by 

(/>„(a;) = V2sm{{n + l)7:x), (27) 

and we set A = Inm. If we again define a„ = 0„(x)dx 
we obtain the set of parameters presented in Table (III). 



TABLE III: Mode dependent vibrational frequencies a;„ and 
integration constants an for the first four streching modes 
with even index of a grapliene nano-ribbon. 



n 









1 


0.90 


2 


3 


0.30 


4 


5 


0.18 


6 


7 


0.13 



As we can see, the a„ decline more rapidly than for car- 
bon nanotubes while the frequencies increase at a slower 
rate. Both these facts would make a graphene nano- 
ribbon inferior to a carbon nanotube if it hadn't been 
for the fact that the amplitude range within the linear 
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FIG. 5: Saturation amplitudes Ao, A2, A4 and Ae for the 
respective modes uo, U2, U4 and Me, as a result of computer 
simulations aimed at selective self-excitation of the second 
harmonic of a graphene nano-ribbon of width w = Inm and 
length L — Ifim. The coupling parameter /3 was in the range 
-9 to -13 which, assuming a stationary current of 20 fiA, cor- 
responds to a magnetic field strength in the interval 6.6-9.6 T. 
The other parameters were Q — 10'', ujl = 1.732, ujb, — 2.000. 
Initial conditions were in all variables. 
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FIG. 6: Saturation amplitudes Ao, A2, A4 and Ae for the 
respective modes ito, U2, «4 and Me, as a result of computer 
simulations aimed at selective self-excitation of the fourth har- 
monic of a graphene nano-ribbon of width Milium and length 
L — l/im. The coupling parameter (5 was in the range -14 to - 
18 which, assuming a stationary current of 20 fj,A, corresponds 
to a magnetic field strength in the interval 10.3-13.2 T. The 
other parameters were Q = 4 ■ 10'', lol = 3.873, lor = 2.000. 
Initial conditions were in all variables. 



regime is greater by almost a factor ten already when 
comparing a graphene nano-ribbon of length 1/xm to a 
nanotube of the same length and with radius 1 nm. This 
can be further extended for graphene by simply making 
it longer. Numerical simulations were performed, now 
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with the couphng parameter 



/3 



HVq tl 
ARq Xp ' 



(28) 



where A is the area of the graphene sheet and to = 
L^y p/Tq. One might think that a small area is bene- 
ficial, but on the other hand a larger sheet has a higher 
current carrying capacity so these two factors are likely 
to cancel each other out. In Figure (5) we show the sat- 
uration amplitudes resulting from a number of computer 
simulations aimed at exciting the second mode. One ob- 
vious difference compared with carbon nanotubes is that 
there is less of a difference between the saturation ampli- 
tude of the excited mode and that of the other modes. 
For graphene they differ typically by a factor 10 whereas 
for nanotubes it is almost by a factor 100. The same con- 
clusion can be drawn from the results presented in Figure 
(6), where the fourth mode is excited. 



harmonic of flexural vibrations of graphene nano-ribbons 
and carbon nanotubes with quality factors of 10* and 
upwards. This is accomplished by the means of a con- 
stant voltage source, an inductor, a capacitor, a gate elec- 
trode and a constant magnetic field. An advantage with 
graphene is that it allows a higher amplitude of oscillation 
before non-linear forces spoils the scheme, but, computer 
simulations indicate that for carbon nanotubes the ratio 
between the saturation amplitude of the excited mode 
and that of the other modes is substantially larger, thus 
allowing a more distinct selective excitation. 
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